> 



The Gravitational Lens System B1422+231: Dark Matter, 
Superluminal Expansion and the Hubble Constant 

David W. Hogg & R. D. Blandford 
130-33 Caltech, Pasadena, CA 91125. 

'rn : 

0\ • submitted 1993 July 23, revised 1993 December 1 

> : 

O , Abstract 

A gravitational lens model of the radio quasar B1422+231 is presented which can account 

■ for the image arrangement and approximately for the relative magnifications. The locations 
\ of the principal lensing mass and a more distant secondary mass concentration were predicted 

and subsequently luminous galaxies were found at these locations. This argues against the 
existence of substantial numbers of "dark" galaxies. The model suggests that if the compact 
. radio source is intrinsically superluminal then the observed component motions may be as large 

\ as ~ 100c with image B moving in the opposite direction to images A and C. The prospects for 

■ a measuring the Hubble constant from a model incorporating lens galaxy locations, compact 
. radio source expansion speeds and radio time delays, if and when these are measured, are 

^ ' briefly assessed. 
^ ■ 

Qh| Key words: gravitational lensing — quasars: individual: B1422+231 — quasars: general — dark 

Q ■ matter — distance scale. 

L- ' 1 Introduction 

•I— I 

^ Gravitational lensing provides a direct probe of the mass distributions in astronomical objects. It is 
' therefore important to find accurate models of lens potentials. Indeed, the existence of a compelling 
lens model can provide corroboration of multiple imaging. In addition, an accurate, verifiable model 
opens the possibility of measuring the Hubble constant if the source proves to be variable. 

The B1422+231 system was discovered by Patnaik et al. (1992) as part of a survey of flat- 
spectrum radio sources. It has a redshift oi z = 3.62, and consists of three bright components (A, 
B and C) and one dim component (D) all within 1.3 arcseconds. There is no evidence, as yet, of 
extended emission. Although D is too dim for accurate spectral and polarization measurements. A, 
B and C have similar radio spectra and fractional polarization. The location of D is consistent with 
the hypothesis that it is the fourth image of a gravitational lens system. Lawrence et al. (1992) 
observed B1422+231 in the infrared and found that all four components are in similar positions and 
have roughly the same flux ratios as in the radio. For these reasons, both sets of authors conclude 
that B1422+231 is a gravitational lens system. 

Simple gravitational-lens models of B1422+231 predicted a mass concentration inside the circle 
of images, plus a mass concentration to the southeast (Hogg & Blandford 1993). Recently several 
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galaxies were discovered in or near the B1422+231 system. Yee & Ellingson (1993) have found a 
galaxy (Gl) inside the circle of images, and Larkin et al. (1993) have found two more galaxies (G2, 
G3) to the southeast of the quasar. 

In this paper, we describe a simple lens model which provides an adequate description of existing 
observational data and predicts the velocity dispersions or masses of the lensing objects, the relative 
directions and approximate relative expansion speeds of the compact radio components (if they turn 
out to be intrinsically superluminal), and the relative time delays between intensity variations in 
the images. We find that using the nearby galaxy locations in our model significantly improves 
our fits, suggesting that the recently discovered galaxies are the main contributors to the lensing 
potential. Future observations (e.g. with VLBl) should improve our ability to make accurate models, 
and thereby improve the remaining predictions. 

2 Observations 

The source B1422+231 was observed by Patnaik et al. (1992) with the VLA at 8.4 GHz, and with 
the MERLIN array at 5 GHz (Table 1). The total flux density is roughly 0.5 Jy at 5 GHz, with 
~ 3% polarization. In the optical, the combined source has the spectrum of a 16.5°^ luminous, red 
quasar, with absolute magnitude My — —29.5 + 51og/i (where h is the Hubble constant in units of 
100 kms~^Mpc~^, and = !)• Observations by Lawrence et al. (1992) in the 2.0-2.4/i K band 
made with the Caltech 58 x 62 InSb array camera at the Cassegrain focus of the Hale Telescope 
are consistent with the radio results. 

The primary galaxy Gl was observed by Yee & EUingson (1993) with the CFHT. When the four 
known components of the B1422+231 system were subtracted, the galaxy Gl was found inside the 
circle of images. G2 and G3 were observed by Larkin ct al. (1993) on a 2/j image taken by the Keck 
Telescope. These two galaxies are very bright; G2 and G3 have K-hand magnitudes of 16.2 and 
15.7 respectively. The positions of Gl, G2, and G3 are given in Table 2. 

A prehminary measurement of the redshift of Gl of z = 0.64 was found by Hammer et al. (1993). 

3 Lens model 

3.1 Gravitational lens theory 

The lens equation relates image positions to source positions. The potential of a gravitational lens 
can be projected onto a two-dimensional plane (the lens plane) orthogonal to the direction of light 
propagation and scaled so that the lens equation is 

/3(0)=0-a(0) = 0-V0V'(0) (1) 

where /3 (a two-dimensional vector) is the angular position of the source (on the source plane), is 
that of an image point (on the image plane), a is the reduced deflection, is the two-dimensional 
gradient operator with respect to 0, and ip{0) is the deflection potential (e.g. Schneider, Ehlers & 
Falco 1992; Blandford & Narayan 1992). 

The inverse of the Jacobian matrix of the mapping /3(0) from the image plane to the source 
plane is the magnification tensor [/x(0)], 

(2) 
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The (scalar) magnification /i of an image relative to the source is the determinant of the magnifi- 
cation tensor. 

The deflection potential ip is related to the "surface" mass density (mass per unit solid angle) 
E by 

VV^^-, where ^c = ^^, (3) 

is the critical density and where D^, Dg, and Dds are the angular diameter distances from observer 

to lens ("deflector"), observer to source, and lens to source. 

There is a time delay At associated with an image at from a source at /3 given by 



At = To 



, where = {1 + z^)^ (4) 



where z^, is the redshift of the lensing potential. By Fermat's principle the images are located at 
stationary values of At{6) (e.g. Schneider et al. 1992). 



3.2 Fitting the data 

We use a general procedure to fit a gravitational lens model to the data (cf. Kayser 1990; Kochanek 
1991a). In the case of B1422-I-231, our data comprise no = 12 observations {0^}, namely the 
positions and fluxes of the four images. We suppose that these observations are normally distributed 
with variances af which we estimate from the results of Patnaik et al. (1992). We then consider lens 
models characterized by um model parameters denoted {Mj}. The number of degrees of freedom 
is 1/ — no — Um (cf. Kochanek 1991b). The model parameters include lens parameters plus three 
more parameters which specify the location and intensity of the source on the source plane. 
We seek a best-fit solution by minimizing a least squares figure of merit 

^^f^mm^ (5) 

with respect to variation of {M,}, where the 0[{{Mj}) are the values for the Oj that we derive from 
the model (e.g. Press et al. 1992). 

Roughly speaking, a fit is good if /y is on the order of unity. Adding additional parameters 
to a model will almost always allow a reduction in x^, but it is only significant if it reduces x^/^- 
This rule of thumb allows us to test the significance of adding additional parameters. 

Prior to performing the fit, it was necessary to explore model space using the time delay surface. 
We identified all the extrema and then minimized a modified obtained by deriving a source 
location and magnification for each image (by backwards ray-tracing), given a model, and then 
using the magnification tensors to convert the source displacements to image displacements for a 
given source location. This is a fast procedure for searching large volumes of parameter space. 

Given a particular model which minimizes x^ we form the Hessian matrix, 

numerically, by varying the model parameters. The inverse of matrix Djji is the covariance matrix 
Cjji which provides standard deviations for the parameters. 
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The goal of this model-fitting is to predict future observables, which we denote by {Pk}- In the 

present case these include the relative magnification tensors of the three bright images and the time 
delays for variations of all four images. We compute the variances in the future observables in the 
usual manner: 



3.3 Lens potential models 

The history of this project is relevant to understanding the results. We first found that a lens 
model based upon a single elliptical galaxy is insufficient to account for the large magnification 
ratio between images A, B, C and image D. The simplest successful lens model, which we call 
the initial model, comprises two singular isothermal spheres with a common redshift, each with 
potential ip — b\0 — 0c\, where the parameter b — Aira'^Df^s/c^Ds, the critical radius, measures the 
ID velocity dispersion a, and Oc locates the centre of the potential (e.g. Kochanek 1991a; Schneider 
et al. 1992). (In practice, having a second potential breaks the circular symmetry and avoids the 
structural instability of its imaging properties.) The initial model has um = 9 parameters, so u = 3. 

On the basis of this model (and simple variants) we predicted that there should be a primary 
lensing galaxy located about 1/3 of the way from image D to image B, and that there should be a 
larger mass located about 7" southeast of image B with equivalent velocity dispersion ~ 460 kms~^; 
presumably a small group (Hogg & Blandford 1993). We then learned that Yee & Ellingson (1993) 
found the primary galaxy at the predicted location. Hammer et al. (1993) made a preliminary 
measurement of the Gl redshift, and Larkin et al. (1993) discovered the two galaxies G2 and G3 to 
the southeast near the predicted location. 

We therefore explored refined models in which we fixed the locations and redshifts of the three 
galaxies Gl, G2, and G3 (we assume that G2 and G3 are at the same redshift as Gl), and adjusted 
their mass distributions. We refer to the best model found this way as the refined model. In 
this case, the two distant galaxies G2 and G3 are modelled as point masses, each with potential 
ip — b'^la\6 — 6c\, where again 6^ locates the centre of the potential, and 6^ = Mg/nHc measures the 
mass Mg of the galaxy (e.g. Kochanek 1991a; Schneider et al. 1992). Treating a distant galaxy as a 
point mass is valid if its mass distribution is roughly circular (on the sky) and none of it overlaps 
the circle of images. The refined model we settle on has hm = 6, i/ = 6. 



4 Results 

4.1 Lens model 

The initial model is the simplest capable of approximately reproducing the observations. The centre 
of the primary potential is located inside the circle of images, and the secondary potential is south- 
east of D. The best-fit parameters together with the errors derived on the basis of the model are 
given in Table 3. In Table 4 we give the computed magnification tensors and time delays for images 
A, C and D relative to image B together. The largest deficiency of this model is its failure to 
reproduce adequately the relative magnification of images A and B. 

We experimented with several sets of adjustable parameters for the refined model. The simplest 
set consists of only three parameters: a velocity dispersion a for Gl, and masses M2 and M3 for G2 
and G3. In this case Gl is modelled as a singular isothermal sphere, and G2 and G3 are treated 
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as point masses. This set of parameters can be expanded by adding a shear to the Gl potential, 
changing the radial dependence of the Gl potential, or by fine-tuning the centre of the Gl potential 
(always remaining consistent with the observations of Ycc & Ellingson (1993). 

We found that introducing the simple, three-parameter refined model reduced x^/t^ from ~ 60 
(for the initial model) to ~ 16. We found that adding parameters to the refined model, such as 
an extra shear, a modified Gl radial dependence, or a fine-tuned Gl location, did not reduce 'X^ jv 
any further. We therefore settled on the very simple three-parameter refined model. The best-fit 
parameters for the refined model are given in Table 3, computed magnification tensors and time 
delays are given in Table 4. Again, the largest deficiency of the model is in reproducing the A/B 
magnification ratio. 

The refined model has a jv which is still quite large (~ 16), but it must be remembered 

that the value depends on the uncertainties in the fiuxes, which we have taken from the radio 
observations alone. The descrepancies between the radio and infrared observations suggest that the 
uncertainties in the fiuxes may be higher than we have assumed. Either a ~ 20% reduction in the 
radio flux of A, or a factor of 2-3 increase in the assumed uncertainties in the A, B, and C flux 
measurements, would be necessary to bring our best-flt jv values to unity. In addition, the lens 
model is very simple and depends upon the assumptions that the galaxies share the same redshift, 
and that there are no significant additional perturbations along the line of sight. 

Changing the variation of mass with radius changes the time delay between the three bright 
images and D. For this reason we believe that the B-D time delay is more uncertain than Table 4 
suggests, though more simulations are necessary to verify this assertion. 

(The ambitious reader may notice that our model produces two spurious images near the centres 
of G2 and G3. These images are only artifacts of our treatment of G2 and G3 as point masses, an 
approximation valid only far from the centres of G2 and G3. In fact, the actual mass distributions 
of G2 and G3 are spread out so that they do not exceed the critical density necessary to produce 
additional images of the quasar.) 

4.2 Lens properties 

In the refined model, the equivalent velocity dispersion of Gl is 210 kms~^ and the mass within 
the critical radius is 9.3 x 10^°/i~^ M©, consistent with a single normal galaxy. Galaxies G2 and G3 

have masses 4.1 x VS^'^hr^ Mq and 4.7 x 10^^/i~^ Mq respectively. These numbers were calculated 
assuming i^o = 1 ^"^(^ -^d = 0.64. If the redshift is significantly smaller than 0.64, the masses 
of G2 and G3 need not be as large to produce the same effect in the model (Browne, personal 
communication) . 

These results are consistent with the predictions of the initial model. This provides an argument 
against the existence of a large density of "dark" galaxies as gravitational lenses ought to sample 
all the mass in the universe fairly. 

4.3 VLBI observations 

The four images of B1422+231 have a combined fiux density of ~ 0.5 Jy. It should therefore be 
mappable by very long baseline interferomctry (VLBl). A VLBI map that resolves the bright com- 
ponents will allow for more detailed models because magnification tensors, not just magnifications, 
may then be included in the fit. 
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If the source has a simple core-jet structure with superluminally moving features, then it should 

be possible to test the hypothesis that the core is stationary. Furthermore, the strong derived 
tangential (i.e. along the A-B-C arc) magnification in our models almost guarantees that the ob- 
served motion in the three bright images will also be tangential (with B's motion reversed). The 
absolute expansion speeds depend upon the unknown intrinsic value but it could exceed ~ 100c or 
~ 9 masyr"^. 

4.4 Time delay 

If additional information can be obtained from radio observations, then it will remove some of the 
uncertainty in the shape of the time delay surface and allow a moderately accurate prediction of 
the relative time delay of images A and C with respect to image B, now anticipated to be several 
hours. (That image B's variation should follow that of images A and C and precede that of image 
D is mandated by the topology of the arrival time surface.) It is possible that delays this short can 
be measured if the source is an intraday variable (e.g. Quirrenbach et al. 1991). In this case, it will 
be possible to fix the multiplicative constant Tq in the time delay and make a measurement of the 
Hubble constant. It may be easier to measure the time delay between images D and B. Measurement 
of relative magnification tensors will also help as they will remove some of the uncertainty in the 
shape of the time delay surface. 

Even if either of these approaches furnishes a measurement of the Hubble constant, there will 
still be a residual uncertainty associated with the specific choice of world model and the possibility 
that a large screen covers the source. For example, changing Qq to ~ 0, while retaining = 0.64, 
increases Tq by 25 percent. Introducing a cosmological constant while retaining fiatness for Qq = 0.1 
(conforming with model C of Carroll, Press & Turner 1992) increases Tq by 12 percent. 

More speculatively, Gott, Park & Lee (1989) used the observations of Q2016-I-112 to exclude 
the possibility of an antipode at large redshift. This argument is much stronger in the case of 
B1422+231 because we can exhibit a successful model and the source is known to be at a higher 
redshift. 

5 Conclusions 

In this paper, we have exhibited a simple gravitational lens model for the radio-loud quasar 
B1422-I-231 despite the unprecedented, large relative magnification of ~ 50. That this is possi- 
ble, and that the quality of the fit improved after the three galaxies were discovered, supports 
Patnaik et al.'s (1992) and Lawrence et al.'s (1992) claim that this source is multiply imaged. The 
primary lensing galaxy Gl is located within the circle of images. The galaxies G2 and G3 act 
perturbatively on the image arrangement, displacing image D closer to Gl (and demagnifying its 
fiux in the process) and having the opposite effect on the other three images. B1422-I-231 provides 
a prime example of magnification bias. At both radio and visual frequencies, the source has un- 
commonly large fiux. The combined magnification on the basis of the refined model is = 29. 
The source quasar is therefore only 20™ and the a posteriori probability of its having been multiply 
imaged is then less striking as there is a much larger parent population. 

As the source is very bright at radio wavelengths it should be possible to obtain high dynamic 
range maps of images A, B, and C. We predict that all three images are elongated along the A-B-C 
arc, and that the parity of image B should be opposite to the parities of images A and C. We 
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also predict large proper motions for the components along the A-B-C arc. Detailed VLBI maps 
that allow measurement of the relative magnification tensors may also lead to more detailed and 
accurate lens models. 

Our final quantitative prediction is the time delay between image variations. If B1422+231 turns 
out to be a strong intraday variable, then it may be possible to measure this for the three bright 
images. (Success in this endeavor would also provide a very clean proof that intraday variability 
is intrinsic to compact radio sources and is not imprinted by refractive interstellar scintillation.) 
Relatively large flux changes in the source will be necessary to measure the predicted ~ 1 month 
delay in the variation of the D image. Time alone will tell if the four images of B1422+231 are 
formed by a simple enough potential and if the source has sufficient intrinsic variability to furnish 
a useful measurement of the Hubble constant. 
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Note added in proof 

Recent radio observations by Browne (personal communication) confirm the A/B fiux ratio which 
our model is unable to fit. M. Remy et al. (1993) have observed the B1422+231 system in several 
optical bands. Their observations are consistent with the morphology found in the radio and IR, and 
they also find galaxies G2 and G3. Because our position for Gl was measured from a preliminary 
sketch provided by Yee & Ellingson, the position we used for Gl differs somewhat from the position 
they have published. The difference is not significant. Kormann, Schneider & Bartelmann (preprint) 
have shown that the masses of G2 and G3 can be reduced if Gl is given an ellipticity. 
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Tables 



Table 1. Image positions. 

image 9i 02 \fJ,^'^^\ 



A -0.39 ±0.02 0.32 ±0.02 0.98 ± 0.02 

B 0.00 0.00 1.00 

C 0.33 ±0.02 -0.75 ±0.02 0.52 ± 0.02 

D -0.94 ±0.03 -0.81 ±0.03 0.020 ±0.005 



Relative positions and 5 GHz radio fluxes of the components of B 1422+231 (Patnaik et al. 
1992). The coordinates 6i, 62, advance in the directions of decreasing RA and increasing 
declination respectively and are measured in arcseconds. 
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Table 2. Galaxy positions. 



galaxy 


Oi 


02 


Gl 


-0.70 


-0.59 


G2 


-9.0 


-5.2 


G3 


-3.6 


-7.3 



Positions of the observed lensing galaxies Gl (Yee & Ellingson 1993), G2, and G3 (Larkin 
et al. 1993), in the same coordinate system as the data in Table 1. 
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Table 3. Model parameters. 



initial model 


refined model 


6(1) 


0.62 ±0.02 




0.75 ±0.01 


^cl 


-0.67 ±0.03 




5.0 ±0.7 


nil) 


-0.60 ±0.02 


6(G3) 


1.7 ±0.3 


6(2) 


3.02 ± 0.02 






(9(2) 


-6.6 ±0.4 






nm 

^c2 


-5.0 ±0.3 







Parameters for the initial and refined models with standard computed 



11 



Table 4. Future observables. 



observable initial model 



refined model 



3.2 3.5 




2.89 3.18 


-0.81 -1.1 




-0.88 -1.23 


" 0.27 0.32 " 




" 0.34 0.37 " 


4.0 2.9 




3.57 2.62 


' 0.65 0.53 " 




" 0.64 0.53 ' 


0.38 0.37 




0.33 0.32 





-0.7 


-0.77 




-0.5 


-0.42 




0.04 


0.03 




-0.2/1-1 


-0.16/1-1 ± 0.04 




-5/1-1 


-1.0/i-i±0.1 




+18/1-1 


+29.9/1-1 ± 0.8 



The computed magnification tensors [/i*^*^^] (and magnifications fi^'^^^) relative to imag 
and time delays At^^^^ in days after B, adopting the fiducial value Tq = 62h~^ d (arcsec)" 
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Figures 

Figure 1. Time delay surface for the refined model. The images are located at A, B, C, D 
and elongated along the same directions as the contours. Contours are separated by about 
l.lh~^ days. 
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